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A measurement of neutrinoless double beta decay in one isotope does not allow to determine 
the underlying physics mechanism. We discuss the discrimination of mechanisms for neutrinoless 
double beta decay by comparing ratios of half life measurements for different isotopes. Six prominent 
examples for specific new physics contributions to neutrinoless double beta decay are analyzed. We 
find that the change in corresponding ratios of half lives varies from 60% for supersymmetric models 
up to a factor of 5-20 for extra-dimensional and left-right-symmetric mechanisms. 
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{N| ' An uncontroversial detection of neutrinoless double beta (OV/3/3) decay fH, 0, 0, 0] will be a discovery of uttermost 
significance. Most importantly, it will prove lepton number to be broken in Nature, and neutrinos to be Majorana 
particles |5|. On the other hand, it will immediately generate another puzzle: what is the mechanism that triggers 
• the decay? The most prominently discussed mechanism for neutrinoless double beta decay is the exchange of light 
l/"") ' Majorana neutrinos. But other mechanisms, like the exchange of SUSY superpartners with R-parity violating or 
conserving couplings, leptoquarks, right-handed W-bosons or Kaluza-Klcin excitations, among others, have been 
discussed in the literature as well. Possibilities to disentangle at least some of the possible mechanisms include the 
analysis of angular correlations between the emitted electrons [3, @] or a comparative study of OV/3/3 and 0v(3 + with 
ly— ^ ' electron capture (EC) decay [7[. Another possibility seems to be the study of double beta decay to excited + states 
Q. Unfortunately, the search for 0vf3 + /EC decay is complicated due to small rates and the experimental challenge 
to observe the produced X-rays or Auger electrons, and most double beta experiments of the next generation are not 
sensitive to electron tracks or transitions to excited states. 

Without identification of the underlying mechanism, an experimental evidence for neutrinoless double beta decay 
. will only provide ambiguous information about the concrete physics underlying the decay. For example, no information 
about the neutrino mass can be obtained from a measurement of the neutrinoless double beta decay half life. 

In general, contributions to neutrinoless double beta decay can be categorized as either long-range or short-range 
interactions. In the first case, the diagram involves two vertices which are pointlike at the Fermi scale, and the 
Oh. exchange of a light neutrino in between, and is described by an effective Lagrangian of the type [t| 

C = ~h (jv-A J v-A,p + ^2 e NpjNpJNP^J , (1) 
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where the sum runs over all Lorentz invariant combinations of hadronic and leptonic Lorentz currents of defined 
helicity, Jnp,v-a — uOjd and Jnp,v~a — eOjV, respectively. Here Oj t j denotes the corresponding transition 
operator. The effective coupling strengths in new physics contributions are denoted as ejvp throughout. For short- 
ranged contributions, on the other hand, the interactions are described by a single vertex being pointlike at the Fermi 
scale. The decay rate therefore results from first order perturbation theory, and is described by the Lagrangian [Io| 

G 2 

£ = —^~ m p 1 ^ ^npJnpJnpj'np- (2) 

Here m p denotes the proton mass and the sum runs over all Lorentz invariant combinations of hadronic, Jnp — uOjd, 
and leptonic, j' NP = eOje c , currents of defined chirality. 

The combination involving two vertices of the first term in JT]) leads to the usual neutrinoless double beta decay 
half life formula for the mass mechanism, 

[T^p 1 = ((m v )/m e ) 2 G 01 \M m »\ 2 , (3) 

where (m„) is the effective neutrino mass in which the contributions of individual neutrino mass eigenstates are 
weighted by mixing matrix elements squared, (m„) = | ^ U^m^. The combination of the first term in (JTJ with any 
of the latter terms as well as the short-range Lagrangian @ leads to the expression 

[T^]-^el P G NP \M NP ?. (4) 
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Here, A4 m " and Ai NP are the nuclear matrix elements for the mass mechanism and alternative new physics contribu- 
tions, and Goi and Gnp denote the corresponding phase space integrals from the list given in We have assumed, 
that one mechanism dominates the double beta decay rate, and we do not consider interference between different 
mechanisms. Calculational details and results for the relevant matrix elements involved have been given elsewhere 
@> E3] ) and numerical results for all common double beta emitter isotopes will be published soon 11 1 . 



In the present context, we will concentrate on the observation that the combinations of lcptonic and hadronic 
currents specific to different mechanisms result in different nuclear matrix elements. This fact taken alone is not 
of much help in order to disentangle the different mechanisms, since e.g. a smaller nuclear matrix element for the 
mass mechanism as compared to any alternative new physics mechanism can be compensated by a larger value for the 
neutrino mass, at least within the constraints implied by other observations such as Tritium beta decay and cosmology. 
However, under the assumption that one mechanism dominates in triggering the decay, the new physics parameter 
{ra v ) or cnp drops out in the ratio of experimentally determined half lives for two different emitter isotopes, 

T 1/2 ( A X) _ \MC e Ge)\ 2 GC 6 Ge) ... 



T 1/2 ( 76 Ge) \M( A X)\ 2 G( A X) ' 

Consequently, half life ratios depend on the mechanism of double beta decay, but not on the new physics parameter, 
and thus can be compared with the theoretical prediction for different mechanisms. Moreover, the error in the 
isotope nuclear matrix element ratio can be reduced compared to the theoretical error in one matrix element, due to 
cancellations of systematic effects. 

In the following we study several prominent examples of specific alternative new physics contributions by calculating 
the corresponding ratios of half lives 

rpNP (A V\ 

v NP,A Y s _ J V2 I ^> (fa 

11 ( X) - T$V6Ge)' (6) 

where we concentrate on a comparison with 76 Ge as it constitutes the best tested isotope to date. We choose the 
following mechanisms for a detailed discussion: 

• SUSY-accompanied neutrinoless double beta decay: 7£ SUSYacc 
This mechanism has been first discussed in |12| . The effective Lagrangian for the dominant contribution assumes 



the form 



£ -j GpU*i ^sUSYacc 

AV2 



(!7 i (l + 7s )e c )(il(l + 75M) + 5(^' i,/ (l+7s)e c )(iI<T^(l + 7 5 )d) . (7) 



and results from integrating out a heavy d-squark of the fc-th generation with i?-parity violating couplings 
A' llfe and A' lfcl , and exchanging a light neutrino of the i-th generation between the nucleons. The new physics 
parameter is given by 



where 9k parametrizes the left-right sfermion mixing of the mass eigenstates d\ and di. 

• Gluino exchange mechanism in R-parity violating SUSY: JZ SVSY ~^ 
In this short-range contribution discussed in [13 , [l^ | , integrating out u- and d-squarks and a gluino leads to the 
effective Lagrangian 

C D ^ m ;V (W+75)d)(«(l + 75)d)-^(T^ (9) 



with 



- 2ira s Xf n m p 



mi 



(10) 



3 





82 Se 


100 Mo 


i28 Te 


130 Te 


136 Xe 
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Ref. 




0.26 


0.11 


3.26 


0.18 


0.77 


0.02 


this paper 


TjSUSYacc 


0.28 


0.11 


3.22 


0.17 


0.53 


0.02 


this paper 




0.28 


0.10 


3.16 


0.17 


0.53 


0.01 


M 




0.29 


0.13 


2.96 


0.20 


0.54 


0.02 
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0.14 


0.13 


18.40 


0.13 


0.67 


0.01 


[17] 


TZ KK (10 GeV- 1 ) 


0.24 


0.08 


3.26 


0.19 


3.31 


0.08 


[15] 


TZ KK (0.1 GeV" 1 ) 


0.26 


0.11 


3.26 


0.18 


0.78 


0.02 


[15] 



TABLE I: Ratios 1Z( A X) of half lives for various important double beta decay emitter isotopes, normalized to the half-life of 
76 Ge. For the exchange of Kaluza-Klein excitations in extra dimensional theories the brane shift parameter and bulk radius do 
not factorize, and are chosen to be a = 10 GeV" 1 , 0.1 GeV" 1 and R = (1/300) eV" 1 . 



• Right-handed currents: n LR - m and n LR ^ xx 

Integrating out right-handed 14^-bosons occurring in left-right symmetric models can lead to two types of new 
contributions with right-handed leptonic currents [H, 

C D S (77i7m(1 + 75)gC) (P^ 1 ~ l5)d) + A( " 7 " (1 + l5)d) ) ' 
where the new physics parameters are given by ry and A. 

• Kaluza-Klein neutrino exchange in extra-dimensional models: 1Z KK 

In extra-dimensional theories, the double beta observable is given by a sum over contributions from all Kaluza- 
Klein excitations with masses m(„), weighted with the mass dependent matrix element Af m "(m( n )) |15j : 



e 



KK _ 1 ^tt2 



U 2 en m {n) {M m «{m [n) )-M m »). (12) 



In this case the effective coupling constant e KK depends on the nuclear matrix element Ai m " (m(„)), and therefore 
the particle physics does not decouple from the nuclear physics. This is because the masses of the Kaluza-Klein 
excitations vary from values much smaller than the nuclear Fermi momentum pf to values much larger than 
Pf, while the TO(„)-dependence of A^ m " (mr n )) changes around pp- Therefore the Kaluza-Klein spectrum has to 
be fixed by choosing specific values for the brane shift parameter a and the radius of the extra dimension 7?. In 
the limit of a — * or R — ► 0, TZ KK approaches TZ m ". 

The matrix elements for the mass mechanism and for the SUSY accompanied neutrino exchange have been calculated 
in the pn- QRPA approach of @, Ell, in the latter case for the first time. For the other mechanisms, existing numerical 
values obtained with the same nuclear structure model have been adopted from the literature. The values for the 
phase space integral factors Goi, Gnp have been calculated in Numerical values for K NP ( A X) are given in Table 
1, and Fig. 1 displays the relative change expected from various new physics contributions, compared to the mass 
mechanism. An application of the procedure to any other alternative new physics contribution by using the matrix 
elements listed in [11| is straightforward. 

All isotope ratios have been normalized to the half life of the most extensively studied nucleus 76 Ge. Moreover, 
while at present no experiment using a 128 Te source has been proposed, we included this isotope since it provides a 
particularly powerful discriminator and thus may encourage future experimental efforts to study this nucleus. 

The two supersymmetric contributions show similar deviations, which are rather small for all isotopes. It is obvious 
that these mechanisms are most effectively discriminated from the mass mechanism by comparing the half life ratios 
between 82 Se and 136 Xe which vary by 60%. In left-right symmetric models, strong deviations can be found for the AA 
combinations, while deviations for the rjrj combination are rather small. A comparison of half life ratios between 100 Mo 
and 136 Xe yields a variation of 70 % for the r\r\ contribution with right-handed hadronic currents, while a comparison 
of measurements in 128 Te and 150 Nd will provide a powerful discriminator with a variation of more than a factor of 
20 for the AA contribution with left-handed hadronic currents. Similarly in extra-dimensional neutrino models with a 
large brane shift parameter, large deviations can be found for 136 Xe and 150 Nd, and the half life ratios for 150 Nd and 
100 Mo vary by more than a factor of 5. Some caution is necessary when referring to the half life ratio of the heavily 
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FIG. 1: Relative deviations of half life ratios 1Z NP ( A X), normalized to the half-life of 76 Ge, compared to the ratio in the mass 
mechanism ft m "( A X). 



deformed 150 Nd, which is ignored in most QRPA calculations (compare the discussion in [18]). Finally it should be 
stressed that not necessarily two positive results are needed - already the comparison of one half life measurement 
and one upper bound in another isotope could provide non-trivial information on the double beta mechanism. 

Since the theoretical errors of the nuclear matrix clement calculation dominate the experimental errors, it is difficult 
to determine the confidence level with which either mechanism can be excluded to generate the observed double beta 
evidence. If, for example, a statistical distribution of matrix element values is assumed, a relative variation of 60% 
in 1Z NP ( A X) with respect to lZ m "( A X) is significant only if the corresponding nuclear matrix elements would be 
known with an accuracy of 15%, which seems to be unrealistic, if only one pair of isotopes is being analyzed. Indeed, 
estimates of errors in nuclear matrix elements vary from a factor 3-5, when the spread of published values is used as 
a measure, to only 30%, according to an assessment of uncertainties inherent in QRPA [TOf . 

However, the significance of the comparison of two isotopes will increase if a whole set of measurements in different 
isotopes resembles the expected pattern. Moreover, one would expect that systematical effects, like an overestimation 
of the nuclear matrix elements due to a too small value for the particle-particle interaction g pp in the pn-QRPA 
approach, a different value for the axial-vector coupling gA, the inclusion of higher-order terms or a different model- 
space would influence calculations for the different isotopes in a similar way, and thereby cancel in the half life ratios 
discussed. This expectation is confirmed by the comparison of the results of different QRPA codes in and of 
QRPA and shell model codes in 3] . Finally it has been pointed out in [13] that the half life ratios ([6|) can also be 
used to single out the correct nuclear structure model. In this case the correct combination of mechanism and nuclear 
structure code can be determined by the best fit of the theoretical half life ratios to half life measurements in various 
nuclei. Thus the results presented in this letter should be complemented and checked with alternative codes for the 
nuclear matrix element calculation. Moreover, other mechanisms, including pion exchange 2l| . may be dominating 
in some of the models discussed, and should be discussed as well. 

In summary, we discussed how different mechanisms of neutrinoless double beta decay would manifest themselves 
in half life ratios involving different isotopes. We thus conclude that complementary measurements in different 
isotopes would be strongly encouraged. At present, next-generation experiment proposals exist for 76 Ge (GERDA, 
MAJORANA, GEM, GeH 4 ), 82 Se (Super-NEMO, DCBA, SeF 6 ), 100 Mo (MOON), 130 Te (CUORE), 136 Xe (EXO, 
XMASS, Xe), as well as for the isotopes 48 Ca, 116 Cd and 160 Gd not discussed in this letter (CANDLES, COBRA 
and GSO) (for recent overviews of the experimental status see [22| ). An experimental study of this kind should 
be complemented by neutrino mass searches in Tritium beta decay experiments and cosmology, as well as searches 
for effects of the alternative new physics source of lepton number violation in other processes, such as lepton flavor 
violating decays . 



After this paper had been submitted for publication, the paper 24J appeared, which comes to similar conclusions 
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and estimates the number of required measurements and their precision needed. 
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